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ABSTRACT 

Protons produced in electromagnetic showers formed by the reverse-electron flux are 
usually the largest component of the time-averaged polar-cap open magnetic-flux 
line current in neutron stars with positive corotational charge density. Although the 
electric-field boundary conditions in the corotating frame are time-independent, in- 
stabilities on both medium and short time-scales cause the current to alternate be- 
tween states in which either protons or positrons and ions form the major component. 
These properties are briefly discussed in relation to nulling and microstructure in 
radio pulsars, pair production in an outer gap, and neutron stars with high surface 
temperatures. 
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1 INTRODUCTION 

Polar-cap models of pulsar radio and X-ray emission all as- 
sume electron-positron pair creation by particles acceler- 
ated along the narrow bundle of open magnetic flux lines 
that intersect the light cylinder. Understanding of the radio 
emission process remains limited. In a review of the sub- 
ject, Melrose (1995) made the interesting observation that 
its characteristics are broadly unchanged over more than 
five orders of magnitude variation in magnetic dipole mo- 
ment. He also noted that whilst instabilities in the electron- 
positron plasma received much attention, there had been 
relatively less interest in the origin of the plasma. But con- 
siderable progress with the latter problem has been made by 
Muslimov & Tsygan (1992) who recognized that the Lense- 
Thirring effect, the general-relativistic dragging of inertial 
frames, has a profound influence on the acceleration electric 
field, which is the component Em locally parallel with B, 
present in the frame of reference rotating with the neutron 
star. The electric field is conservative in the rotating frame 
and is given in terms of the charge density p by, 



V • E = 4tt (p - Po ) , 
where, in Euclidean space, 



po 



4-kc 



V • ((O X r) X B) 



(1) 



(2) 



is the density at which force-free corotation of the magne- 
tosphere exists (Goldreich & Julian 1969). The most sig- 
nificant way in which a general-relativistic treatment (see 
Muslimov & Harding 1997; equations 37 and 38) modifies 
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equation (2) is that the rotation angular velocity SI seen by 
a distant observer is replaced by SI — Slur, where Qlt is 
the Lense-Thirring angular velocity. In the non-relativistic 
equation (1), both p and po have the same radial dependence 
because the angle between SI and B does not change much 
over distances no more than an order of magnitude greater 
than the polar cap radius. But p and po have different phys- 
ical origins, p being determined principally by continuity, so 
that the presence of the Lense-Thirring angular velocity 



2G 



(3) 



where L is the angular momentum of the star, in the general- 
relativistic expression for po produces a non-trivial radial 
dependence of the quantity p — po which is the source of the 
acceleration field. 

In the corotating frame, the potential boundary condi- 
tion for equation (1) is $ = on the neutron star surface 
and on the surface separating open from closed magnetic 
flux lines, with E = — V<&. This is supplemented by the 
condition Em = on the neutron star surface if the work 
function there is not large enough to affect the free flow of 
charge. Thus there are three possible cases: (i) SI • B > 0, 
po < 0, E|| = with electron acceleration; (ii) SI ■ B < 0, 
po > 0, Ey = 0, giving ion and positron acceleration or 
(iii), SI ■ B < with Ey ^ 0. The latter case is that consid- 
ered originally by Ruderman & Sutherland (1975) in which 
there is now renewed interest (see Gil et al 2008, and work 
cited therein) owing to the possibility that some neutron- 
star polar-cap magnetic fields may exceed the critical field 
B c = m 2 c 3 /eTi — 4.41 x 10 13 G and may even be one or 
more orders of magnitude larger than the dipole moment in- 
ferred from spin-down indicates. Recent calculations of ionic 
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work functions (Medin & Lai 2006) , including magnetic flux 
densities B > 10 14 G, do not unequivocally exclude the 
neutron-star surface boundary condition Ey 7^ 0. 

Following the work of Muslimov & Tsygan, case (i) ap- 
pears consistent with electron acceleration and pair produc- 
tion occurring over a very wide range of surface magnetic 
flux densities. Ion acceleration occurs in case (ii) with the 
possibility of positron production by secondary processes. 
However, the time-independent boundary conditions $ = 
and E|| = applying in cases (i) and (ii) are very restric- 
tive and in terms of electromagnetic considerations alone 
there is no obvious reason why either the current density or 
the nature of the electron-positron plasma within the open 
magnetic flux-line region should change so as to produce 
phenomena such as pulse nulling over times of the order of 
10 1 — 10 5 s. (For example, failure of the condition $ = on 
the surface separating open from closed magnetic flux lines 
over such long time intervals appears unlikely.) In this con- 
nection, it is worth noting that in PSR B1931+24, nulls last 
for intervals so long that it has been possible to measure 
the spin-down rate separately during the on and off-states 
of emission (Kramer et al 2006). Spin-down in the on-state 
is approximately twice as fast as in the off-states. In each 
of cases (i) and (ii), the self-consistent current density in 
the open magnetic flux-line region (see Muslimov & Hard- 
ing 1997) does not deviate by so large a factor from the 
Goldreich- Julian value of poc. Thus some change in the na- 
ture of the plasma is indicated. Also there is the fact that the 
radio emission from many pulsars exhibits short-term vari- 
ability in the form of both microstructure and subpulses. 
The microstructure may be caused by temporal variation of 
the accelerated plasma or may be a consequence of angular 
beaming, but there appears to be as yet no conscensus (see, 
for example, Smirnova et al 1994 and Lange et al 1998). It 
is not obvious how any of these phenomena arise from elec- 
tron or positron acceleration subject to time- independent 
boundary conditions. For these reasons, the present paper 
addresses the question of whether or not it is plausible that 
case (iii) actually exists in certain radio pulsars and, specifi- 
cally, to an investigation of nuclear processes in the electro- 
magnetic showers produced by the reverse flow of electrons 
toward the polar cap. In the course of this, we also find 
that case (ii) is actually less simple than might be supposed 
and that, even in the presence of time-independent bound- 
ary conditions, natural instabilities exist that cause medium 
and short time-scale variability in the nature of the acceler- 
ated plasma. 

Nuclear processes were considered previously at mag- 
netic flux densities B < B c in a less than satisfactory at- 
tempt to estimate backward photon production by the elec- 
tromagnetic showers that are formed at the neutron-star sur- 
face (Jones 1979). But the present work is concerned with 
B ~ B c because only in this region is there a possibility 
that case (iii) could exist. Also the problem here is much 
more simple: that of finding the rate of photo-production 
of protons by decay of the nuclear giant-dipole resonance 
(GDR), for which experimental data has greatly improved. 
The mean electron-shower energy required to create a pro- 
ton is a function of B and of the nuclear charge Z, but at 
B > B c is found to be no more than of the order of 2 — 5 
GeV. Diffusion to the neutron-star surface occurs with a 
short, but important, time delay and the process leads di- 



rectly to an upper limit for a steady-state surface blackbody 
temperature in cases (ii) and (iii). But further analysis of the 
short-time instability in case (iii) shows that T res , the polar- 
cap temperature in the absence of a reverse-electron flux, is 
actually the temperature which should be compared with 
the critical temperature T c inferred from calculated ionic 
work functions. The details of shower development and pro- 
ton formation and diffusion leading to these conclusions are 
given in Sections 2 and 3. The application to pulsar polar 
caps is described in Section 4. 

Proton production by the reverse flux of electrons is 
prolific even at B <C B c (see Section 2.5) and in Section 
4.3 it is shown that it generally leads to medium time-scale 
instability in cases (ii) and (iii) producing similar variabilty 
in electron-positron pair formation. Instability with a short 
time-scale is also likely to be a property of cases (ii) and (iii) 
but unfortunately, is not possible to predict with complete 
confidence because it depends on details of the proton dif- 
fusion to the polar-cap surface that are not known with cer- 
tainty. However, it is of some interest in view of the observed 
characteristics of individual radio pulses in some pulsars. 

The principal results of this paper are the proton forma- 
tion rates at B ~ B c and the instabilities found in the fluxes 
of accelerated nuclei and protons and in pair formation. Sec- 
tion 5 considers briefly the possibility that these may be con- 
nected with pulse nulling and microstructure, pair formation 
in any outer-magnetosphic gap, and the radio-observability 
of young neutron stars with high whole-surface tempera- 
tures. 



2 ELECTROMAGNETIC SHOWER 
DEVELOPMENT 

The later stages in the development of an electromagnetic 
shower whose primary particle of energy Ea ~ 10 2 - 10 3 
GeV moves parallel with B are determined principally by 
Coulomb pair production, Compton scattering, and by mag- 
netic bremsstrahlung activated by multiple Coulomb scat- 
tering. Formation of the giant-dipole state is by far the most 
important nuclear interaction and occurs predominantly in 
the very late stages of shower development, within the pho- 
ton energy band of 15 — 30 MeV. Higher-multipole interac- 
tions have much smaller cross-sections and appear only at 
greater energies where the photon track length within the 
shower is smaller. The principal processes will be described 
in some detail in Sections 2.1 - 2.4, but in Section 2.4, other 
processes that are present in showers at B > B c are dis- 
cussed briefly. 

The natural unit of length for shower development at 
B <C B c is the radiation length defined in terms of the 
Bethe-Heitler formulae (Bethe 1934) for the bremsstrahlung 
and pair production cross-sections. This convention is 
adopted here for convenience, though with cross-sections 
modified by a screening wavenumber defined by the mat- 
ter density at the neutron-star surface. For 10 < B12 < 10 3 , 
where B12 is the magnetic flux density in units of 10 12 G, wc 
fit the ion number densities found by Medin & Lai (2006) 
by the expression 

N = 2.7 x W 25 Z^? 7 Bl 2 2 cm' 3 . (4) 
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The radiation length is then, 



2.4 x W~ 2 Zn 



(5) 



in which the logarithm, whose argument is derived from 
the infinite linear-chain ion separations given by Medin & 
Lai, replaces the ln(f 83Z -1//3 ) of the Bethe-Heitler formulae. 
This is extremely small; l r = 5.2 x 1CT 5 cm at B = 10 14 G 
and Z = 26. The Bethe-Heitler mean free path for Coulomb 
pair production is 9l r /7- 

The electron or positron dispersion relation in a field of 
magnetic flux density B is 

, p., 1/2 

E=[l+p\+2n—J . (6) 

Unless otherwise stated, all electron or photon energies and 
momenta will be expressed in units of the electron rest en- 
ergy and of mc. The classical momentum components py 
and p± are parallel with and perpendicular to B, and n is 
the Landau quantum number. The conserved quantities in 
electron-photon interactions are the total energy and the 
parallel momentum component, but not p± except in the 
correspondence principle limit of n S> 1. States of n ~ 1 
are important at B > B c , and are best regarded as of 
one-dimensional motion, given by pii , with effective mass 
yjl + 2nB/B c . It is unfortunate that there appear to be no 
published calculations of the Coulomb bremsstrahlung and 
pair creation cross sections that are valid for B ~ B c . 



2.1 Magnetic bremsstrahlung 

Multiple Coulomb scattering of a shower electron or positron 
is a sequence of ascending (or descending) Landau transi- 
tions from its creation state n±. For B > B c , the clas- 
sical (n ^> 1) energy loss-rate expression for magnetic 
bremsstrahlung (see the reviews of Erber 1966 and of Hard- 
ing & Lai 2006) is not valid, particularly in the later stages 
of shower development where the dominant process has two 
stages: Coulomb excitation of the n — — > 1 transition 
followed by cyclotron emission. In order to see that this is 
so, and in the absence of an accurate expression, we are 
obliged to adopt the correspondence principle for an order 
of magnitude estima te and t herefore assume that a momen- 
tum transfer q± ~ \J2B/ B c from the nuclear Coulomb field 
is necessary for an n — — > 1 transition. The mean free 
path, derived from the Rutherford scattering formula for q± 
and from the number density of equation (4), is then, 



A = 5.1 x 10~ J Bi 2 hi 



(62Zi 



/2 R -l/2 
26 13 12 



(7) 



equivalent to ~ 10 2 Coulomb excitations per radiation 
length at 10 14 G. The relativistic quantum-cyclotron emis- 
sion rate can be approximated by t^} ~ 1.0 x 10 _3 a;s in the 
py = frame of reference, where u>b is the classical cyclotron 
angular frequency. This expression is almost stationary with 
respect to n and B/B c in the vicinity of 10 14 G (see Fig. 3 
of Harding & Lai 2006). The condition that the cyclotron 
emission rate exceeds the Coulomb excitation rate is that, 



7o < 



0.7S?2 8 Z 2 - 6 13 , 



(8) 



in which 70 is the Lorentz factor for transformation from the 
Pll = frame to that of the rotating star. Thus the statement 



that n = ^ n = 1 transitions are the dominant processes 
is valid in the very late stages of shower development at 
10 14 G and we can conclude that in this region electrons or 
positrons lose their energy over distances one or two orders 
of magnitude smaller than l r . 

From the kinematic conservation laws, the maximum 
photon transverse momentum (that is, _!_ to B) in cyclotron 
emission for ann-+n' transition is, 



k± = \ / 1 + 2n 



B_ 

~B C 



1 +2n' 



B_ 

K 



(9) 



Because the low-energy part of the electron-positron spec- 
trum is the source of most GDR-band photons, we shall 
show, in Section 2.3, that this limit on photon transverse 
momentum is of crucial importance in estimating their total 
track length in the shower. But earlier in the shower, for 
values of 70 about two orders of magnitude larger such that 
7ocr ce > l r , Coulomb bremsstrahlung is the more important 
energy-loss process. 

2.2 Coulomb bremsstrahlung 

The cross-section for Coulomb bremsstrahlung at B > B c 
must differ from the Bethe-Heitler formula, but no result 
appears to have been published at the present time. For the 
emission of a photon of momentum k at an angle 6 with B, 
by an electron or positron of initial momentum py with a 
change of Landau quantum number n — ► n' , the required 
momentum transfer components are, in the limit py >• 1, 

B 1 

q« = fc(l - cos#) + (n - n)— h 

Bc P\\ 



k ( , B\ 

^{ 1 + 2n BZ)> 



q± 




-—, for n n 



(10) 



For transitions with n' 7^ n, the threshold longitudinal 
momentum transfer is larger than its zero-field value of 
gy = fc/2p 2 and it is therefore possible that the true cross- 
section is reduced in comparison with the Bethe-Heitler 
value. But this merely increases the electron track length 
in the earlier stages of shower development, with which we 
are not concerned, and does not affect the photon total track 
length in the GDR band. 



2.3 Coulomb and magnetic pair creation 

A photon with momentum k at an angle 9 with B can 
convert to an electron-positron pair with Landau quantum 
numbers n±. Only the longitudinal component of momen- 
tum transfer from the nuclear Coulomb field appears in the 
kinematic conservation laws and the threshold is, 

qn = k{ l- coa e)-J-(l + 2n + §-)- 

which is valid in the limit py ± S> 1 and has two regions 
of interest; gy = and gy / 0. It is again unfortunate 
that there appears to be no published cross-section for the 
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Coulomb process at finite q» , and in finding the total photon 
track length in this case, we shall assume the Bethe-Heitler 
mean free path, multiplied by an unknown correction fac- 
tor 7] p > 1. However, q\\ — is possible at a finite 6, and 
the associated transverse momentum component k± is the 
threshold for single-photon magnetic pair production. The 
individual thresholds that are important here are, 

k± = 2, n+ = n_ = 0, 



k± = f + yt + 2— , n±=0, n T = l. (12) 

Transitions to the Landau state n_ = n + = are not al- 
lowed for photons polarized so that the electric vector is 
parallel with k x B (Semionova & Leahy 2001, see also Hard- 
ing & Lai 2006), but are allowed for photons polarized with 
electric vector perpendicular to k x B. The classical expres- 
sion for the conversion rate, which is valid in the n± ^> 1 
limit, is a function of k±B/B c (see Erber 1966; Harding & 
Lai 2006), and is not applicable here at fields B > B c with 
n± ~ 1. But the rate R je e = 1.2 x 10 6 cm -1 that it pre- 
dicts for 30 MeV photons at 10 14 G, immediately above the 
higher of the two thresholds given by equations (12), is so 
large that, whatever its inadequacies, it must be presumed 
that magnetic conversion occurs within a very small frac- 
tion of a radiation length. Thus photons with k± above the 
higher threshold do not contribute significantly to the total 
photon track length or to the photoproduction of protons 
through GDR formation. Photons in the allowed polariza- 
tion state can convert above the lower threshold but the rate 
is reduced by a factor of the order of (k±/k) 2 w 10~ 3 and 
is approximately an order of magnitude smaller than that 
for Coulomb conversion. A further small contribution to the 
magnetic conversion rate must exist because the correct elec- 
tron and positron propagators are those for passage through 
condensed matter, not free space. Therefore, Coulomb inter- 
action introduces a small component of n± — 0, n T = 1 into 
the n+ — n- = state. But this is relatively unimportant 
and has not been studied in detail here. 

Comparison of equations (9) and (12) shows that all 
photons emitted by cyclotron decay of Landau states n to 
the ground state n — do not exceed the thresholds for sub- 
sequent single-photon magnetic pair production provided, 

\l 1 + 2n w c < 3 ' 

\l 1 + 2n K < 2 + ^+¥ (13) 

where the two inequalities refer to the two thresholds of 
equation (12). The angular distribution of the cyclotron de- 
cay photon in the pu = frame of reference considered in 
Section 2.1 is approximately isotropic (see Fig. 3 of Latal 
1986), so that because equation (9) gives the maximum pos- 
sible value of k±, photons with k± below the maximum may 
well fall below the thresholds (12) even if the above inequali- 
ties are not satisfied. It is interesting to note the first inequal- 
ity fails for the minimum Landau quantum number n = 1 
for all B > 4B C whilst the second is satisfied for all B. 



2.4 Compton scattering and other processes 

Compton scattering is the most important of the processes of 
second order in the fine structure constant that are not con- 
sidered in previous Sections. It can significantly modify the 
very late stages of shower development at B > B c . Extensive 
cross-section calculations have been performed by Gonthier 
et al (2000) for a wide interval of photon energy. For the 
purposes of estimating its significance, the mean free path 
equivalent to the Thomson cross-section ot for the electron 
density given by equation (4) with B = 10 14 G and Z = 26 
is At = 8.5 x 10~ 6 cm, which is less than a radiation length. 
The total cross-sections (see Fig. 2 of Gonthier et al) are 
negligible at photon energies k > 10 3 , but are significant in 
the GDR band. Scattered photons with k± exceeding the 
higher threshold (12) undergo prompt magnetic conversion. 
Gonthier et al have shown that the exact cross-section in 
this region is well approximated by the non-magnetic Klein- 
Nishina formula, a very convenient result which enables us 
to establish that scattered photons with k± below the mag- 
netic pair thresholds comprise only a very small, and for 
present purposes, negligible fraction of the total Compton 
cross section. Thus the cross-sections that effectively convert 
photons to electron-positron pairs can be obtained directly 
from Fig. 2 of Gonthier et al. They are slowly varying func- 
tions of B/B c . We shall assume a photon energy of 21 MeV 
for GDR formation in nuclei of atomic number 10 ^ Z ^ 26 
(see the review of Hayward 1963). At this energy, the Comp- 
ton cross sections are ac = 3.6 x 10~ 2 <tt for B = B c and 
4.8 x 10~ 2 (tt for B = 10B C . At higher energies, the cross 
section decreases, being approximately oc fe _1 . 

The other second-order processes which should be con- 
sidered briefly involve the blackbody radiation field, whose 
photons are typically of energy k bb ~ 10~ 3 . Inverse Compton 
scattering can give photons exceeding the pair production 
thresholds (12), but only for electron momenta greater than 
k ~ 10 4 . Similarly, two-photon pair production is possible, 
but with a threshold k — k^ 1 and a negligible transition 
rate because the blackbody photon number density is many 
orders of magnitude smaller than the ion densities given by 
equation (4). Direct production of electron-positron pairs 
(the trident process) is of third order and has no signifi- 
cant effect on shower development. Photon splitting, also 
of third order, has a small transition rate, saturating at 
R ps « 0Ak _1 k± cm -1 at B > B c (see the review of Harding 
& Lai 2006). 

The remaining process that does change shower devel- 
opment is the Landau-Pomeranchuk-Migdal effect for which 
we refer to the extensive review given by Klein (1999) , also to 
the work of Hansen et al (2004) for its laboratory verification 
at electron energies of several hundred GeV. In condensed 
matter at zero field, multiple Coulomb scattering removes 
the coherence that the Bethe-Heitler formulae assume over 
distances of the order of g^ 1 , where is the very small 
longitudinal momentum transfer from the nucleus, and so 
reduces the cross sections, particularly for bremsstrahlung 
production of low-energy photons. This is almost certain to 
be a feature of the showers considered here owing to the high 
density of matter at the neutron-star surface. But it is not 
possible to make even qualitative estimates of the effect be- 
cause the nature of multiple Coulomb scattering is changed 
at B > B c (see Section 2.1) and there is also a strong one- 
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dimensional ordering of ions parallel with B arising from the 
formation of bound moleculsr chains (see Medin & Lai 2006) 
which will be of significance for questions of coherence of the 
Coulomb bremsstrahlung and pair creation amplitudes. It is 
almost certain that in the early stage of shower develop- 
ment considered here, the LPM effect leads to an inward 
displacement of the GDR photon band region, though with- 
out significant change in its total track length. We shall see 
in Section 3 that the formation of protons is also not sig- 
nificantly changed but that their time-scale for diffusion to 
the surface is modified in a non-trivial way as described in 
Section 4.3. 



2.5 Photon track length in a shower 

The distribution of photon track length is found from elec- 
tromagnetic shower theory. On dimensional grounds, the 
track length per unit interval of photon energy k such that 
E c <C k <C Eo in a shower of primary electron energy Eo 
must be of the form, 



G ( k > = —gr 



(14) 



in which y is a dimensionless parameter and E c is the crit- 
ical energy at which electron or positron ionization and 
bremsstrahlung energy-loss rates are equal. In the zero-field 
case, where there is more certainty in the fundamental pro- 
cesses of shower development, it is possible to obtain reli- 
able values of y. But at B > B c , we are obliged to adopt a 
more elementary and analytical approach, usually referred 
to as Approximation A, which includes only the most signif- 
icant processes which are bremsstrahlung and pair creation 
(Landau &: Rumer 1938). But even at this level, there are 
difficulties because the relative importance of Coulomb and 
magnetic bremsstrahlung changes as the shower develops. 
The later stages of shower development with which we are 
concerned here are almost completely determined by mag- 
netic bremsstrahlung, as described in Section 2.1. Therefore, 
a modified Approximation A, in which throughout the devel- 
opment of the shower, Coulomb bremsstrahlung is replaced 
by magnetic bremsstrahlung, has been used to calculate an 
approximate value of y. Pair creation is assumed to be de- 
termined by the Bethe-Heitler formula with the unknown 
correction factor r\ v retained (Section 2.3). To obtain the 
total photon track length, only the form of the photon spec- 
trum is needed and initially, we assume the number density 
tx fc~ 2 / 3 per unit interval of k which is valid in the classical 
magnetic bremsstrahlung limit of n ^> 1. Following closely 
the procedures of Landau & Rumer (1938) and Nordheim 
& Hebb (1939), we find that G(k) is given by the inverse 
Mellin transform, 



(15) 



in which 
D 



2 8 1 



3s+2 3s 
7 / 1\ f 1 £ s - 1 



C) 2/3 



(16) 



with ( chosen so that all poles lie to the left of the contour. 
To obtain the term linear in Eo we retain only the contribu- 



tion from the s = 1 pole and find y — 0.608%,. Apart from 
the unknown factor, this differs little from the zero-field Ap- 
proximation A in which, for an approximate photon number 
density oc fc _1 , the value is y — 0.47. 

But as shown in Section 2.1, the magnetic 
bremsstrahlung classical limit is not reached in the 
later stages of shower development. Instead, the dominant 
processes are Landau transitions n = 0^n=lin which 
the angular distribution of the cyclotron decay photons in 
the jjii = frame of reference is approximately isotropic. 
Thus a Lorentz transformation to the frame of the rotating 
star produces a photon number density distribution which 
is a constant per unit interval of k. In this case, the 
denominator function D(s) in equation (15) is replaced by 

_ . , 2 8 1 8 1 7 
Dc(s) = 1 s, (17) 



s+1 3 s + 2 3s +3 9 

giving the result y — 0.871%. It is worth noting that the Ap- 
proximation A results for total photon track length are in- 
dependent of the absolute rates of bremsstrahlung emission. 
Even the effects of quite substantial changes in its spectral 
shape are relatively small, being limited to small increases 
in y for the harder spectra. 

Two processes not included in Approximation A, Comp- 
ton scattering and energy loss by ionization, can affect the 
late stages of shower development. The principal effect of 
Compton scattering is a reduced effective mean free path for 
pair creation by photons in the GDR band and therefore, a 
decreased total photon track length available for GDR for- 
mation. This is important and can be allowed for by the 
replacement, 



Vp -> T~~ . x c = ^Nl r Zac, (18) 

1 + T) V XC 7 

in equation (15), with ac(k) evaluated at the GDR forma- 
tion energy k g d — 41. 

It is known through comparison of Approximations A 
and B that in the zero-field case, electron and positron ion- 
ization energy loss produces some reduction in y because 
GDR-band energies are of the same order of magnitude 
as the critical energy E c . But the total electron-positron 
track length in the later stages of shower development at 
B > B c is much reduced because energy loss by magnetic 
bremsstrahlung predominates and occurs over lengths one 
or two orders of magnitude smaller than l r (see Section 
2.1). Thus the relative importance of ionization energy loss 
is much reduced near the GDR band and we attempt no 
correction for it. 

On this basis, we shall adopt the track-length estimate 
given by equation (14) in which yl r is found from equations 
(4), (5), (15), (17) and (18), with the parameter r\ v as an 
unknown. 



3 PROTON PRODUCTION 

The mean number of protons formed in a shower per unit 
interval of primary electron energy is 



W p = {NL lp )G(k gd )/E , 



(19) 



in which 7 7P is the energy-integrated cross section for giant 
dipole resonance formation and decay into channels contain- 
ing one or more protons, and the square brackets denote an 
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average over the depth distribution of GDR-band photons 
in the shower. This latter function is unknown in detail, 
owing to the gradual transition from Coulomb to magnetic 
bremsstrahlung as the dominant process in shower develop- 
ment and to the existence both of the unknown parameter 
r]p and the LPM effect, but it is certainly very small near 
the surface and increases to a maximum at a depth which is 
likely to be of the order of 10i r - Therefore, we are obliged to 
evaluate the ion number density N and I 1P at a mean atomic 
number Z. The integrated cross-section is expressed as the 
electric-dipole sum-rule value (Hayward 1963) multiplied by 
a factor x p , 

I-y P = 117.4 (x v Z) (l - mc mb (20) 

where A is the nuclear mass number. Measured cross- 
sections for light elements 10 < Z < 26 provide evidence for 
a value x p w 0.5 (Hayward 1963; Kerkhove et al 1985; As- 
safiri & Thompson 1986). But in the circumstances of shower 
development below the neutron-star surface, the problem of 
estimating the effective value of x p is more complicated. 

The protons formed in GDR decay lose energy by ion- 
ization, are quickly thermalized, and diffuse to the surface 
with negligible probability of secondary nuclear interaction. 
However, this is not so in the case of GDR-decay neutrons 
which lose energy and diffuse to the surface by a sequence of 
nuclear elastic scatterings in the course of which capture in a 
transition A — > A + 1 may occur with a significant probabil- 
ity which is, however, hard to estimate. Consequently, nuclei 
in the shower region may be neutron-rich. There is some evi- 
dence (Fultz et al 1974; O'Keefe et al 1987) that x p decreases 
with increasing nuclear asymmetry A — 2Z, possibly owing to 
the associated divergence between neutron and proton sepa- 
ration energies. A simple model of the effect of this increase 
in nuclear asymmetry is as follows. The energy-integrated 
cross-section for neutrons is given by equation (20) with the 
substitution p — > n and Xp ~\~ Xfl — Q> w 1. Suppose that neu- 
trons diffuse to the surface with survival probability e n and 
that the captured fraction 1 — e n increases nuclear asym- 
metry. This produces a bias toward GDR decay by neutron 
emission. To take account of this, let x n — Hx p , in which H 
is a function of A — 2Z. The steady-state condition for equal 
neutron and proton loss rates is then = Ht n Xp, 

from which we find that e n determines the asymmetry func- 
tion H and x p = ae„(l + e„) _1 . In the limit of small e n , 
the nuclear asymmetry is, of course, limited by /3-instability, 
the relevant time-scale being that in which the corotational 
current density poc removes the equivalent of one radiation 
length of matter, t H = Nl r AeP/B ~ 10 2 s for B = 10 14 G 
and rotation period P = 1 s. 

Values of W p found from equation (19) are given in the 
final column of Table 1 for intervals of Z and of B that 
might permit the Ey =fc boundary condition at the polar 
cap. They assume the parameters r\ v = 1 and x v = 0.5. To 
allow for the effect of neutron capture, the numbers in the 
right-hand column should be multiplied by 2ae n (l + e n ) _1 . 
Multiplication by r) p (l + xc){l + r/pXc)^ 1 allows for values 
of rip > 1. In the latter case, it is worth noting that even in 
the limit 7/ p ^> 1, the existence of the Compton effect, which 
has been extensively studied at B > B c , maintains shower 
development. Proton formation rates in the Table do not 
vary greatly as functions of Z or of B and are extremely 



Table 1. The number of protons created per unit primary 
electron energy, in units of mc 2 , are given in the right-hand 
column for values of the mean nuclear charge Z in the region 
of the shower maximum and of the magnetic flux density B in 
units of the critical field B c = 4.41 X 10 14 G. Also given are 
the radiation length l r in neutron-star surface matter and the 
parameter xc that allows for the Compton effect correction 
given by equation (18). It is assumed here that the correction 
factor applied to the Bethc-Heitler pair creation mean free 
path is rip = 1 and that the strength of the giant dipolc- 
rcsonance proton decay channels is given by x p = 0.5. We refer 
to Section 3 for the correction to bc applied given different 
values of these two parameters. 
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2.0 x 10~ 4 


1.30 
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10 


9.2 x 10~ 5 


2.96 


2.8 x 10~ 4 


18 
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2.0 x 10~ 4 


0.37 


1.5 x 10~ 4 


18 
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7.3 x 10- 5 


0.58 


1.7 x 10~ 4 


18 


10 


2.9 x 10~ 5 


1.11 


2.2 x 10~ 4 


26 


1 


1.1 x 10~ 4 


0.23 


1.0 x 10- 4 


26 


3 


4.0 x 10~ 5 


0.36 


1.2 x 10~ 4 


26 


10 


1.5 x 10~ 5 


0.64 


1.5 x 10~ 4 



high, so that a primary electron of 10 3 GeV creates between 
200 and 500 protons. It is, of course, the case that many 
approximations have been made in Sections 2.1 - 2.4 but 
these do not have a great effect on the late stages of shower 
development in which most of the total track length given 
by equation (14) is generated. It is principally for this reason 
that we believe Table 1 provides a reliable estimate of the 
true rates. 



3.1 Proton diffusion 

Diffusion of protons to the surface is very rapid. The dis- 
tribution of proton formation depth depends on a number 
of factors and is not easily obtained from shower calcula- 
tions such as those described in Section 2.5. The dominant 
process of photon creation changes during shower develop- 
ment from Coulomb to magnetic bremsstrahlung (see Sec- 
tions 2.1 and 2.2). The LPM effect, considerd in Section 2.4, 
displaces shower development inward as does the unknown 
factor r/p > 1, particularly the region of low-energy photon 
formation which is relatively compact. On this basis, it is 
convenient to assume a normalized proton formation depth 
distribution given by 

g p (z) = 8(z - z v ), (21) 

though with the reservation that z p is not well known. Let 
the time interval between creation and arrival at the surface 
be t — t' — Tt p with r dimensionless. For a semi-infinite 
homogeneous medium, the probability per unit interval of r 
that a proton created at depth z < arrives at the surface 
z = is, 
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in which z p = 2>Dt p and D is the diffusion coefficient. Thus 
the probability of arrival at the surface with delay r is, 

f P (r) = f dzg p (z)P(r,z) (23) 

J — OO 

= (^r 3 ) ' exp(-A). 

The solid at the polar-cap surface consists of linear molecu- 
lar chains of ions, parallel with B and normal to the surface, 
probably ordered so as to form a three-dimensional lattice. 
An approximate upper limit for D can be obtained from the 
proton thermal kinetic energy and from the ion separation 

q 1/2 1/2 

a s fs 5 x fO Z 2 g B 12 cm found from the calculations of 
Medin & Lai (2006). It is, 

D = a °\[^* 4 x W ~\IW^ cmV1 > ( 24 ) 

V P™v V S 12 

where /3 _1 = feT and Te is the polar-cap surface tempera- 
ture in units of 10 6 K. This is based on the reasonable as- 
sumption that, within the solid structure, the proton poten- 
tial barriers are at most ~ fcsT. For z p — 10l r and Te — 1, 
the characteristic time for diffusion to the surface is then 
t p w 10~ 5 with the values of l r given in Table 1. 

The condition of the surface at such temperatures is 
that it remains solid. The melting temperatures predicted 
by the one-component Coulomb-lattice condition given by 
Slattery, Doolen & DeWitt (1980) are close to 10 7 K at B = 
B c and are consistent with the melting temperature shown 
in Fig. 1 of the paper by Potekhin et al (2003) . (We shall see 
in Section 4 that the longitudinal thermal conductivity is so 
large that the same statement can be made with confidence 
of matter at the shower maximum z — z p .) Protons reaching 
the surface are accelerated parallel with the magnetic flux 
density B either immediately or, if their flux exceeds the 
corotation value, from the top of the gravitationally-bound 
proton atmosphere that is then formed. 

It is also worth noting that the evaporation of ions at 
the surface is determined purely by the thermal excitation 
of lattice degrees of freedom. Direct interaction with the in- 
ward flux of photons produced by curvature radiation from 
reverse electrons is not significant in this respect. Photo- 
disintegration of nuclei within a small number of lattice 
planes immediately below the surface has negligible prob- 
ability. Photoelectric or bound-free cross-sections for the 
removal of electrons from ions are large, but the momen- 
tum transfer to the nucleus is determined by the Fourier 
transform of the bound electron wave function and it s max- 
imum value is therefore of the order of l^Tha^ 1 U B/B c , 
where ao is the Bohr radius. The equivalent recoil energy 
is 140 A -1 BB~ X eV for an ion of mass number A, which is 
negligible. 

The magnitude of the photoelectric cross section is itself 
a source of reverse electron flux (see Jones 1981). Electrons 
removed from partially-ionized atoms accelerated through 
the polar-cap blackbody radiation field form a significant 
flux at the polar-cap surface. Given the large values of W p , a 
direct consequence is that the time-averages of the electron- 
positron, ion and proton current densities J e , Jz and J p 
on open magnetic flux lines satisfy J p 3> Je,z- The time- 
averaged particle flux above the polar cap consists almost 
entirely of protons. 



3.2 The structure of the proton and ion 
atmosphere 

Ion cohesive energies are by no means negligible at B » B c . 
Those calculated by Medin & Lai (2006) for Z = 26 can 
be expressed approximately as E c — 0.16£>i 2 7 keV in the 
interval 10 < B12 < 100 which is of interest here, and at 
the lower end of this interval are in fair agreement with the 
cohesive energy of 0.92 keV for Z = 26 at B12 = 10 calcu- 
lated by Jones (1985) using a different model for the three- 
dimensional lattice. Thus f3E c S> 1 and, in consequence, 
there is a large density discontinuity at the transition be- 
tween the solid phase and the very thin gravitationally- 
bound atmosphere. At lower magnetic fields, such that 
PE C ~ 1, the discontinuity is less obvious and the atmo- 
sphere more dense. 

The part of this atmosphere that is in local thermo- 
dynamic equilibrium exists within < z < zi, where we 
assume that 21 is the surface of last scattering, roughly de- 
fined here as the altitude beyond which number densities 
are so small that scattering of electrons and ions can be 
neglected. The proton number density in this interval is in 
general many orders of magnitude smaller than those of ions 
and electrons which therefore determine the atmospheric 
structure. Let us assume for the purposes of illustration that 
there is present only one ion species and that it is partially 
ionized with charge Z and mass number A. The standard 
relations between number densities and chemical potentials 
for free Boltzmann particles (allowing for non-zero B in the 
case of electrons) , with the constraint of electrical neutrality 
(strictly charge density po), require that in local thermody- 
namic equilibrium there must be a very small electrostatic 
potential $(2) given by 

e*(z) = (M AJ - m) gz, (25) 

Z + 1 

where g is the local gravitational acceleration. This result is 
not much changed if allowance is made for the presence of 
different ionization states (see Jones 1986). The scale height 
of the atmosphere is very small and, in the general case, 
the structure is of successive layers of increasing Z/A as the 
altitude 2 increases. Equation (25) can be used to obtain 
the scale height of the very small proton component. It is 
negative, showing that protons cannot be in thermodynamic 
equilibrium at 2 < 21 and are accelerated across the surface 
of last scattering. 

The proton number density equivalent to the corota- 
tional charge density is given by N p ev p /c = po and, even for 
polar-cap thermal values of velocity v p , is much smaller than 
ion densities at z\. Thus as particles move out in Landau 
states parallel with B to 2 > 21, the condition of electrical 
neutrality again requires an electrostatic potential that is 
independent of the small proton number and is also given 
precisely by equation (25). The proton potential energy at 
2 > 21 is 

m p (z - z x )g + e*(z) - e*(zi) = 

m p (z-zi)g- -J— (M A Z -m) (z - zi)g, (26) 
Z + 1 

which is evidently negative. The relatively small numbers 
of protons are again accelerated whilst ions, with smaller 
charge-to-mass ratios, are retained by the gravitational field. 
Thus the protons are preferentially selected for the process 
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of inertial acceleration that follows. If the flux of protons at 
z = is greater than is needed for the corotational current 
density, the surplus forms an atmosphere at z > zi. This 
ordering has significant consequences for plasma formation, 
as will be seen in Section 4. 

3.3 Reverse electron flux from photoelectric 
ionization 

Inertial acceleration under the Ey = polar-cap surface 
boundary condition was first described by Michel (1974). 
Owing to current conservation, an excess of positive charge is 
present immediately above the surface in the non-relativistic 
stage of ion acceleration. This is the source of an electric 
field, 

Ei w (**«g*y = 5.0 x 10* (-^;y°^) 1/2 (27) 

expressed in electrostatic units (esu) at z > Zi, which is the 
roughly defined altitude at which ion motion becomes fully 
relativistic, 

^(^CV /2 = 4.1X10> ( ~ AP V /2 cm( 28 ) 
\18tt p Ze J \B 12 Z cosV / 

where P is the rotation period and ip is the angle between 
CI and B. (In obtaining these approximate results, the ion 
charge Z has been assumed constant in the interval < z < 

Zi.) 

The change of ionization from partial (Z) to complete 
(Z) is the result of photoelectric absorption of blackbody 
photons whose energies have been boosted by Lorentz trans- 
formation to the rest frame of the accelerated ion. A sim- 
ple model of the process has been described previously in 
which transitions occur sequentially in order of increasing 
electron separation energy at altitudes for which unit optical 
depth of blackbody radiation field for a particular transition 
is reached. In this way, a compact approximate expression 
for the mean reverse electron energy at the polar-cap surface 
per accelerated unit positive charge can be obtained (Jones 
1981; see equations 15 - 22, and Table 2 of that paper). 
The prefactor in the expression for the electron separation 
energies assumed there was obtained for smaller magnetic 
fields (B\2 ~ 1) than those considered here, but comparison 
with the recent calculations of Medin & Lai (2006) for mul- 
tiply ionized C and Fe ions shows that separation energies at 
B > B c can found with sufficient accuracy for present pur- 
poses by scaling from that expression. Following this proce- 
dure, the mean electron energy per accelerated unit nuclear 
charge is estimated to be 

(^) C_e " 5 - 6 X lQiZ ' lf ' (°)^ 5T 6" 10 ™> (29) 

in which Z(0) is the mean nuclear charge at the surface 
z = 0. These values assume the acceleration field is uniform 
at altitudes small compared with the polar-cap radius, but 
are quite insensitive to variations in its strength because 
photoionization occurs quite promptly when the blackbody 
radiation field is boosted, in the rest frame of the accelerated 
ion, to the electron separation energy. The altitude at which 
this occurs has little effect on the reverse-electron energy 
flux. We assume a value E\\ — Ei — 10 6 esu. Equation (29) 



obviously fails for very light ions, Z < 4 — 5, which are likely 
to be completely ionized in the LTE region at z < z\ . It can 
be combined with an approximate parameterization of the 
values for W p given here in Table 1, 

W v = 7.0 x 10- 6 {Z^- 76 ) sm B^ 2 (me 2 )" 1 , (30) 

where square brackets denote the average nuclear charge in 
the region of the shower maximum, to obtain an expression 
for the number of protons produced per unit positive nuclear 
charge accelerated. It is, 

K = K Z°- 85 (Q) J g v {z)dz 

in which, 

K = 0.24(Z 26 ' 76 ) smJ B 1 2 62 T 6 - 1 '°. (31) 

We note again that this is invalid for very light ions, Z < 
4 — 5. It is also not valid at much smaller fields where K 
tends to an asymptotic value. Here, W p is given by the zero- 
field value of the shower parameter y (see Section 2.5) and 
by placing B12 = 1, so that equation (30) is replaced by 
W p = 3.8 x 10" 5 (Z : f 6 - 76 ) 3m (mc 2 )" 1 . Electronic separation 
energies concerned here approach asymptotically their zero- 
field values at B12 < 1 so that equation (29) should be 
evaluated at B12 = 1 to obtain an estimate of the asymptotic 
value of Ko, 

Ko = Koas = 0.13(^26° ' 76 }smT 6 (32) 

Equations (31) and (32) are an under-estimate of K 
owing to our neglect here of a further source of reverse elec- 
trons that is difficult to quantify, specifically pair creation 
following inverse Compton scattering of blackbody photons 
by the photoelectrons or the consequent conversion of curva- 
ture radiation photons at higher altitudes. The electron and 
photon momenta are so oriented that the necessary electron 
Lorentz factor, 10 4 — 10 5 , is easily reached. For the signif- 
icance of the outward accelerated positrons as a source of 
pairs, we refer to Harding & Muslimov (2002). 



4 APPLICATION TO POLAR CAPS 

Neutron stars with spin direction such that fi • B < and 
polar-cap surface boundary condition Ey = must presum- 
ably exist but the extent to which they are present in the 
observed population of radio pulsars is an unsolved problem. 
The physical existence of case (iii) of Section 1, the boundary 
condition Ey / 0, is less certain in the observed neutron star 
population and is a question to which proton production in 
polar-cap reverse electron flow is very relevant. This section 
attempts to estimate the interval of polar-cap magnetic field 
that would support this condition and to investigate the na- 
ture of the plasma formed. It does not attempt to consider 
the actual radio emission process. 

The distribution of the inferred polar-cap dipole field 
strengths for the observed radio pulsar population contained 
in the ATNF Catalogue (Manchester et al 2005) indicates 
that a large fraction are likely to be too small to support 
case (iii). But it is interesting to note, though with reserva- 
tions about the possibility of unknown selection bias, that 
the nulling radio pulsars listed in Tables 1 and 2 of the pa- 
per by Wang, Manchester & Johnston (2007) have a distri- 
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bution significantly displaced toward larger fields in com- 
parison with the complete ATNF catalogue population. The 
distribution of log 10 B, where B is here the listed polar- 
cap field, for the whole population of 1486 radio pulsars is 
slightly skewed toward low values but has a mean of 11.98 
and a median of 12.06. With neglect of the low-B tail, the 
distribution has a standard deviation of approximately 0.50. 
The 63 pulsars for which Wang et al give a finite nulling 
fraction have a mean of 12.41 and a median of 12.45. The 
expected standard deviation of the mean for a randomly se- 
lected set of 63 pulsars is approximately 0.07 which shows 
that the set of nulling pulsars have field strengths that are 
significantly larger than the general population. 

We shall see that if the polar-cap fields were identical 
with their inferred dipole values, few if any pulsars would be 
able to support the boundary condition Ey 7^ 0, given the 
ion cohesive energies referred to at the end of Section 3.2. 
But in this argument, it is implicit that the dipole origin is 
positioned at the centre of the star. Outward displacement 
of the origin along the dipole axis toward the surface of the 
star would produce considerable amplification of the field at 
one polar cap and diminution at the other without changing 
significantly the dipole field at the light cylinder. Beyond 
this simple case there are, of course, many other possible 
configurations with enhanced polar-cap fields. 

The proton formation rates given in Table 1 are high 
and any atmosphere formed under boundary condition (ii) 
is fractionated in the manner described in Section 3.2. Thus 
the protons undergo inertial acceleration first and, if ex- 
hausted, are followed by ions in order of decreasing Z/A 
if they are present in the atmosphere. Because the time- 
averaged total open magnetic flux-line current density con- 
sists almost entirely of protons and must satisfy the condi- 
tion J < poc it follows that the values of W v lead directly to 
an upper limit for the reverse electron-photon energy flux at 
the polar-cap surface. The temperature limit derived from 
this for a typical value W p mc 2 = 2 x 10~ 4 is, 

T_=6.2xl0 5 (-^2^) 1/4 K. (33) 

Because it has been obtained from the time-averaged energy 
flux, it can in principle be exceeded in a time-variable state. 
But in a time-independent state it can be compared with the 
critical temperature T c above which the ion thermal emis- 
sion rate is high enough to maintain the case (ii) boundary 
condition Ey = 0. In terms of the work function, this is ap- 
proximately ksT c — 0.030-E c . The values of Medin & Lai 
(2006) for Z = 26 and 10 < B 12 < 100 can be expressed as 
E c = 0.16Bi° 2 7 keV and give T c = 5.6 x 10 4 B? 2 7 K. This can 
be compared with T ma x or with T res , the polar-cap black- 
body temperature in the absence of reverse-electron flux. 
(This is approximately the whole-surface blackbody temper- 
ature corrected to the local proper frame.) We can see that, 
regardless of the reverse electron flux, the Ey boundary 
condition can be satisfied only in very old pulsars or in those 
where the polar cap field is of the order of B c . 

4.1 Medium time-scale variability 

The observed field-strength distribution for nulling pulsars 
indicates that it may be worth considering whether or not 
the variability time-scales associated with nulling can be a 



feature of the boundary conditions assumed in case (iii). 
But first, in Section 4.1, we shall consider case (ii). As we 
have seen in Section 3, photodisintegration reduces nuclear 
mass numbers in the region of the shower maximum. Protons 
diffuse to the surface with negligible probability of nuclear 
interaction. A fraction e n of neutrons also escape nuclear 
capture and reach the surface. Thus nuclei fixed in a lattice 
plane move gradually toward the surface with a small ve- 
locity v. (Our reference to a lattice is merely a convenience. 
There will be order parallel with B, but also a considerable 
density of point defects produced by neutron scattering.) 
Instability against the growth of curvature radiation pairs 
requires a minimum Eu and hence, broadly, that the cur- 
rent density J should not exceed a critical value J cr it- We 
consider first case (ii) with current density J > J cr u so that 
there is no instability against pair formation by curvature 
radiation. Solution of equation (1), with time-independent 
boundary conditions and inclusion of inertial acceleration, 
determines J which is constant. But intuition suggests that 
the distribution of nuclear charge with depth is not neces- 
sarily time-independent because the surface charge Z(0, t) is 
determined by the surface state at a previous time t — t sm , 
where t sm is the time interval in which nuclei move from the 
shower maximum to the surface. Let us assume, for simplic- 
ity, that e„ = 1. 

It will be convenient in this Section to describe the mo- 
tion of nuclei toward the surface in terms of a Lagrangian 
depth z which moves with velocity v(z,t) relative to z so 
that z = z at t = 0. The rate of change of nuclear charge 
within an element at this depth containing a fixed number 
of nuclei, NSz where N is given by equation (4) , is 

87 

NSz^ = ~g p (z)SzKoN(0,t')Z 1+ »(0,t')v(0,t'), (34) 

where v = 0.85 from equations (31) and (32). But N oc Z a 
so that the evolution of Z satisfies, 

Z a+1 {-oo) - Z a+1 {0,t) = 
(q + 1) I dt'5(t-t' -t sm )K Z- +a+1 (0,t'). (35) 

J — oc 

This expression contains the approximation g p = S(z — z p ) ~ 
v _1 (0, t')5(t — t' — t S m) which is satisfactory provided g p is 
a sharply-peaked function displaced below the surface as is 
expected owing to the nature of shower development and to 
the LPM effect. It is also unimportant that t sm is to some 
extent a function of Z- values in the interval z p < z < 0. The 
initial value of Z is Z( — 00) and it is assumed constant. A 
time-independent solution of equation (35) exists, given by, 

( z(o) y +a 1 

\Z(-oo)J l + (l+a)K 

But assuming Z( — 00) = 26, and a time-independent charge 
at the surface so small, Z(0) = 5, that equations (29) and 
(31) begin to fail because few accelerated ions have any 
bound electrons, we find a self-consistent value K — 5.5 
for a = —0.7. However, this state is not stable because for 
a charge fluctuation 8Z(0,t), equation (35) becomes a ho- 
mogeneous Volterra equation of the second kind with no 
non-zero square-integrable solution. The charge fluctuation 
satisfies, 

6Z(0,t) = -(v + a + l)K Z v (p,t-t 3m )SZ(0,t-t am ), (36) 
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and for the above value of K, grows indefinitely with alter- 
nating sign. Ultimately, equation (36) fails and the system 
alternates between intervals of high nuclear charge at the 
surface and intervals in which the surface nuclear charge is 
so low that the reverse electron flux is negligible. The basic 
unit of time in which t sm is measured is the time t r i, de- 
fined in Section 3, for the removal of one radiation length of 
matter at the corotational current density. This is certainly 
within the orders of magnitude associated with nulling in- 
tervals. The behaviour in time of a complete polar cap is, of 
course, much more complex because the interval of depth z p 
considered here is extremely small compared with the cap 
radius uq, and the question arises of the extent (if any) of 
correlation between different areas of the cap. 

It is obvious that there must be variability with time- 
scales of the order of t rt in case (iii). Although there is little 
information about the Z-dependence of the ion work func- 
tion E c , we can be confident that it must ultimately decrease 
to a very small value in the limit Z — ► 1. Thus there exists a 
critical value of Z at which the boundary condition Ey / 
fails locally in some region of the polar cap and this must 
ultimately be reached as a consequence of proton emission. 
But the acceleration field is a function of the condition of 
the whole polar cap and the overall level of complexity is 
such that we make no attempt here to proceed further with 
this analysis except to note that time-independent case (iii) 
solutions are not possible. 

4.2 Factors relevant to short time-scale variability 

There are several characteristic times that are important 
factors in examining the possibility of short time-scale vari- 
ability of plasma production in cases (ii) and (iii). The first 
is the proton diffusion time t p derived from z p and D. Then 
the initial e-folding growth time for a curvature-radiation 
electron-positron cascade must be of the order of t ee = 2l a /c, 
where l a is the open- magnetosphere acceleration length in- 
terval concerned. The third parameter is t con d the time-scale 
for thermal diffusion from the shower maximum z v to the 
surface. 

For an estimate of t con d, we require the surface tem- 
perature fluctuation from T to T + ST(0,t) caused by a 
fluctuation SX(z',t') in the thermal power input at depth 
z < below the polar-cap surface at time t' . The correct 
Green function for this particular problem is cumbersome 
(see Carslaw & Jaeger 1959), but the characteristic value of 
t — t' is, for present purposes, satisfactorily approximated 
by t cond — CZp/2\\\. Here Ay is the longitudinal coefficient 
of thermal conductivity and C is the specific heat of sur- 
face matter, for which an upper limit given by the Dulong 
and Petit law is assumed. Transport coefficients at B > B c 
in neutron-star envelopes and crusts have been calculated 
by Potekhin (1999). The longitudinal coefficient of thermal 
conductivity is not a particularly rapidly varying function 
of B and T and, at Z = 26, B = 3B C and T = 10 6 K, 
is approximately 6 x 10 14 erg cm -1 s _1 K _1 . With the ion 
number density given by equation (4) and a shower depth of 
z p = 10l r , the time is t cond f^5x 1CP 10 s, which is short com- 
pared with t ee and with the estimates of t p obtained from 
equations (21) - (24). It is also possible to confirm that, 
for the expected power inputs associated with the derived 
temperature limit given by equation (33), the temperature 



gradient between the shower maximum and the surface is no 
more than of the order of 10 6 K cm , which is negligible. 

Thus in relation to variability on time-scales set by t ee 
or t p , we can assume that the surface temperature fluctu- 
ations are instantaneous functions of the reverse-electron 
power input to the polar-cap surface. This is one of three 
important factors. The second is the mass-number fraction- 
ation of any atmosphere which was described in Section 3.2. 
The final factor is that the accelerated protons produce a 
quite negligible reverse-electron flux, as can be seen by ex- 
amining their electromagnetic interactions with the thermal 
radiation field. Protons accelerated to energies of the order 
of 10 3 GeV are well below the threshold for pair-creation, 
directly or by Compton scattering, and in any case the cross- 
sections concerned are small. However, accelerated ions are 
also below the pair-creation threshold but, as we have seen 
in Section 3.3, those with Z sufficiently high to have bound 
electrons produce a reverse-electron flux by a sequence of 
photo-dissociation reactions with blackbody photons which 
increase the ion charge from Z toward Z (Jones 1981; see 
Medin & Lai 2006 for atomic ionization energies at B « im- 
possible sources of the electron-positron current density 
J e in neutron stars with tt ■ B < need consideration and 
we refer to the recent review of Harding & Lai (2006). The 
growth of curvature radiation cascades requires a large -Ey , 
that is, a partial vacuum above the polar cap (J < Jcrit), 
with the possibility that J crit = 0. A further source is in- 
verse Compton scattering of the blackbody radiation field 
by the reverse electron flux from accelerated ions, followed 
by magnetic conversion of the scattered photons. This was 
mentioned briefly at the end of Section 3.3 and we refer to 
Harding & Muslimov (2002) for further details. 



4.3 Short time-scale variability; E\\ = 

A first investigation of short time-scale variability in case (ii) 
with T m ax > T c is most easily made, as in Section 4.1, with 
the assumption of a time-independent J, with negligible J e . 
The proton current is then, 

J p (t) = K [ dt'f p (t-t')(j-J p (t')), (37) 

J — oo 

where here the parameter A" is a constant given by equations 
(31) or (32). The time t in this Section is dimensionless and 
in units of t p . Equation (37), with J determined by solu- 
tion of equation (1) for the immutable boundary condition 
$ = on the surface separating open from closed magnetic 
flux lines, may not be a completely adequate description of 
short time-scale variability of the system, but is a useful 
starting point. It has the obvious time-independent solution 
J p — KJ I (K + 1) but it can be seen that for any small 
fluctuation in proton current density, equation (37) reduces 
to a homogeneous Volterra equation of the second kind hav- 
ing no square-integrable non-zero solution. For the present 
problem as opposed to that of Section 4.1, we choose to in- 
troduce a short time-scale fluctuation as an inhomogeneous 
term SJ p o{t) centred on t = 0, and rewrite the equation as, 

Jp(t) = KJ + SJ p0 (t) ~K drf p (r)J p (t - r), (38) 
Jo 
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with r = t — t' . We can see that solutions that are not 
square- integrable exist for K > 1. For example, the case 
fpij) = 5(t — to) gives the obvious iterative solution, 

K 



J p (t) 



K + l 



J ± 8J p o(t) — K8J p o(t — to) + 
K 2 8 J p0 (t-2T )-.... 



(39) 



More generally, we can introduce the Fourier transform of 
J p o into the formal iterative solution of equation (38) and 
sum to infinity to obtain, 



duJ i + Kf P H ' (40) 



where, 



drfpir) exp(iwr). 



(41) 



As expected, the trivial case of equation (38) can be recov- 
ered by summation over the residues of those poles in the 
integrand that are in the upper half of the complex-u; plane. 
On the other hand, there exist simple functions, such as 
f p (r) = acxp(— or) with constant a that give no such poles. 
This indicates that the behaviour of / P (r) at small r > 
is crucial, so that for non-trivial functions, such as equation 
(23), there appears to be no alternative to a numerical search 
for zeros in the denominator. Thus for f p given by equation 
(23) we find, for K = 100, poles at w = ±4.85 + 1.88i. For 
smaller K they approach, but do not reach, the real axis; 
for K — 25, the poles are at ui = ±3.4 + O.li. 

Our conclusion is that there is no doubt that the pos- 
sibility of instability exists for large K with exponential 
growth of short time-scale fluctuations until the condition 
< J p < J is no longer satisfied and equation (37) fails. 
But unfortunately, certainty is not possible because we do 
not have sufficient knowledge of the actual shape of the for- 
mation depth distribution g p at depths z v < z < but 
close to the surface. The existence of poles in the upper half 
of the complex-oj plane in equation (40) appears to depend 
crucially on the time interval between the creation of pro- 
tons and their arrival at the surface. This almost certainly 
depends on the LPM effect considered in Section 2.4. 

The behaviour of the system when equation (37) fails 
is easy to see. Suppose, for example that T res > T c , as must 
be the case in many pulsars, so that the boundary condition 
E |l = is always maintained, also that J cri t < J so that 
growth of curvature radiation is not possible. Instability re- 
sulting in excess proton production leads to the formation 
of a proton atmosphere at z > z\, as described in Section 
3.2, and the surplus results in a current density containing 
only the proton component of magnitude J p = J. For the 
interval of time within which this state persists there is no 
reverse-electron flux, no further proton formation, and the 
surface temperature falls promptly to T res in a time no more 
than an order of magnitude greater than t CO nd- At the end 
of this interval, the proton atmosphere is exhausted and J p 
falls abruptly to a value given by the tail of the distribution 
of f p , much below its time-independent value, KJ/(K + 1). 
Consequently, there will be a prompt burst of ion accelera- 
tion and of outward moving positrons produced by the in- 
verse Compton scattering process (see Harding & Muslimov 
2002) unless T res < T c . Owing to the time-delay, this can 
be expected to lead to the formation of a further proton ex- 



cess. Strict periodicity is not to be expected because, as in 
the case of medium time-scale variability, there is no obvious 
reason why the state of a very thin layer, of height z ~ z\ 
should be uniform over the whole polar-cap area. 



4.4 Short time-scale variability; Ey ^ 

In a time-independent state, this boundary condition re- 
quires T max < T c , where T max is the polar-cap blackbody 
temperature given by equation (33). But in a state of short 
time-scale variability, it can be satisfied intermittently pro- 
vided the less restrictive condition T res < T c is valid. We can 
begin by considering a steady state with J > J cr it as in the 
previous section, with polar-cap temperature T > T c . The 
instability leads to a state of excess proton production and 
to a very prompt decrease in temperature to T rea in a time 
probably one or two orders of magnitude larger than t con d 
but still much smaller than t p and to a similarly fast partial 
collapse of the ion LTE atmosphere, so that only protons 
and any \ow-Z ions remain. When these are exhausted, the 
current density J promptly decreases toward zero and if J cr it 
has no finite value (the star cannot support pair creation by 
curvature radiation at any current density) it would appear 
that plasma generation ceases. But a finite J cri t allows pair 
creation when J < Jcrit- Thus J e + Jz increases rapidly and 
a further proton excess is generated. 

The conclusion to be drawn here is that, given short 
time-scale variability, the important polar-cap surface tem- 
perature is T r cs which should be compared with the work- 
function related temperature T c = 5.6 x 10 4 -Bi 2 7 K. Thus for 
T res — 5 x 10 5 K, a surface magnetic flux density -B12 > 23 is 
required for the boundary condition Em / to be satisfied. 



5 CONCLUSIONS 

There seems to be no known reason why neutron stars with 
i~i ■ B < and hence positive polar-cap corotational charge 
density should not be formed at a rate of the same order as 
those with the opposite sign. The purpose of this paper has 
been to investigate some of the features of plasma formation 
at the polar cap to see if these correlate with the properties 
of any subset of the observed population of radio pulsars. 
Specifically, we have attempted to find the formation rate 
for protons in the electromagnetic showers produced by the 
reverse electron flux at the polar cap. These protons dif- 
fuse rapidly to the surface and, having both a negligibly 
small work function and the highest charge to mass ratio, 
are preferentially accelerated upward by the electric field 
component parallel with B. Proton formation is necessarily 
accompanied by evolution of the nuclear charge within the 
very thin surface layer in which the showers occur. 

The proton formation rate was estimated previously 
(Jones 1981) at magnetic flux densities of B12 ~ 1 but the 
present paper has been addressed to the problem at much 
higher fields (B ~ B c ) owing to the realization that such 
values, though infrequent, are not unknown in the distribu- 
tion of inferred dipole-field strengths. Additionally, an im- 
portant subset, those radio pulsars exhibiting the nulling 
phenomenon, have higher than average fields, as noted at 
the beginning of Section 4. Values of the parameter K, de- 
fined as the number of protons formed per unit positive ac- 
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celerated nuclear charge, are given in Table 1 for B ~ B c , 
but its parametrization in equation (31) is not valid at small 
values of B for which its asymptotic value K as is given by 
equation (32). 

Proton formation and the evolution of nuclear charge 
have the consequence that particle fluxes under the case (ii) 
boundary condition defined in Section 1 display much more 
complexity than might have been supposed. In particular, 
there is an instability, given by equation (36), which is re- 
alized in almost any radio pulsar, giving time variability 
typically over intervals of the order of t r \. This is the time 
in which one radiation length of matter flows from the polar 
cap and is within the orders of magnitude associated with 
nulling intervals. Owing to the proton formation-depth dis- 
tribution, fluctuations in the surface nuclear charge grow un- 
til the system alternates between having high values close to 
Z{— oo) and nuclear charges so low that the reverse electron 
flux is negligible. The latter state has no source of electron- 
positron pair creation at the polar cap and may possibly 
be associated with the intervals of null emission. But the 
state of the whole polar cap is almost certainly much more 
complex than a simple transition between high and low sur- 
face nuclear charge. The total depth of the shower forma- 
tion region and of the LTE atmosphere is very many orders 
of magnitude smaller than the polar cap radius «o ~ 10 4 
cm. The question then is why the time- variable state of the 
surface and atmosphere should be in phase over the whole 
polar-cap area. No serious attempt has been made here to 
investigate this problem but there seems to be no reason 
to expect such phase uniformity. Chaotic behaviour appears 
more probable. Similar time-variability and complexity are 
expected in case (iii). 

Instability leading to short time-scale variability is also 
possible. The relevant unit of time is t p , the characteristic 
time for proton diffusion from depth z v to the surface. Here, 
growth of the instability results in alternation between two 
states of which the first is an accelerated proton current 
density J p = J with J e = Jz = 0, no reverse electron en- 
ergy flux and a polar-cap surface temperature equal to T res . 
Abrupt exhaustion of the proton atmosphere then allows 
rapid growth in J e and Jz and the formation of a short pulse 
of electron-positron pair formation. This instabilty appears 
to be common to both cases (ii) and (iii). It is also interest- 
ing because it defines the temperature at which the case (iii) 
boundary condition / must be satisfied. This must be 
T res and not the temperature T max of equation (33) defined 
by the time-averaged maximum reverse-electron energy flux. 
The condition is then T re3 < T c . From our parametrization 
of the work functions found by Medin & Lai (2006), the 
critical temperature is T c = 5.6 x 10 4 Bi 2 7 K, and we find 
that the required polar-cap surface magnetic flux density is 
B12 > 2.3T r 1 e 4 5 in which T re s5 is the polar-cap temperature, 
in units of 10 5 K, in the absence of a reverse-electron en- 
ergy flux. This unit of T is not inappropriate for the whole- 
surface temperatures of the older radio pulsars and the result 
shows that the case (iii) boundary condition can be realized 
at much lower magnetic flux densities than was previously 
thought. At the same time, the similarities between the in- 
stabilities found in cases (ii) and (iii) suggests that, in terms 
of observable phenomena, the distinction between the two 
boundary conditions may not be of great significance. It is 
to be emphasized that, in both cases, the cessation of ion 



flow is not temperature-controlled (by the ionic analogue of 
Richardson's formula) . Rather, it is controlled by excess pro- 
ton production and by the consequent formation of a proton 
atmosphere at z > z\. In case (iii), only the onset of ion flow 
is a result of surface heating. 

An interesting state is that of case (ii) with the pa- 
rameter K < 1 derived from the B\2 <C 1 expression for 
the asymptotic limit given by equation (32). This can be 
realized at polar-cap temperatures greater than 10 6 K, par- 
ticularly if this temperature is supported by a large value of 
Tres- With reference to equation (36), we can see the asso- 
ciated proton and ion current densities are stable and time- 
independent. Monoenergetic ion beams do interact with an 
electron-positron plasma (see Cheng & Ruderman 1980) but 
we anticipate little if any electron-positron pair production 
by inverse Compton scattering at B12 <C 1 and there may be 
no observable consequences of such a system. This may be 
particularly relevant to young neutron stars or to any older 
neutron stars that for other reasons have very high values 

Of Tres- 

A further comment can be made about the reverse elec- 
tron energy flux from any outer gap that exists. It is hard to 
make quantitative estimates of this, but there is every pos- 
sibility that the consequent rate of proton formation would 
completely inactivate polar-cap pair formation as in nulling 
intervals. 

The time- variabilities considered so far have been found 
assuming a time-independent boundary condition $ = on 
the surface separating open from closed magnetic flux lines. 
But we have not addressed the problem of how, or if, this 
boundary condition is maintained. The proton diffusion time 
t p is little more than an order of magnitude greater than the 
typical transit time of order uoc~ and a correct treatment 
of polar-cap phenomena may require not just the solution of 
equation (1) but solution of the full set of Maxwell equations 
with retardation. Thus we do not exclude the existence of 
further instabilities beyond those considered here. This com- 
ment may also be relevant to case (i) in which the corota- 
tional charge density is negative. The considerations of this 
paper have no application to this class of neutron star. The 
electron work function is so small that there is no doubt 
the condition Em = is satisfied at all times on the polar 
cap, also that electrons are freely available to maintain the 
$ = boundary condition on the surface separating open 
from closed magnetic flux lines. 

It would be of interest to know how boundary conditions 
(i), (ii) and (iii) are distributed within the observed popula- 
tion of radio pulsars and neutron stars and whether or not 
the short time-scale temporal variability shown here to be 
a feature of cases (ii) and (iii) contributes to the observed 
microstucture of some radio pulses. Perhaps systematic mea- 
surements of individual pulse structures with 10~ 6 — 10~ 5 s 
resolution would provide some insight. 
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